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SUMMARY 


Laminar film condensation under the simultaneous influence of gas - 
liquid interface shear and "body force (g force) is analyzed over a flat 
plate. Important parameters governing condensation and heat transfer of 
pure vapor are determined. Mixtures of condensable vapor and nonconden- 
sable gas are also analyzed. The conditions under which the body force 
has a significant influence on condensation are determined. 


INTRODUCTION 


Condensation is one of the basic processes which are important in 
the operation of a power plant. Most of the condensation processes which 
take place on the earth depend on the strength of the gravitational field. 
If condensation is to be carried out effectively in a space vehicle, away 
from a strong influence of gravity, another force should be provided 
which would continuously remove the condensate from the condensing 
surfaces. One of the natural existing forces available for the condensate 
removal in forced flows is the shear force at the gas -liquid interface. 

Nusselt, in his classical analysis of film condensation on a vertical 
plate under the influence of the gravitational force (ref. l), neglected 
the shear stress which exists at the gas -liquid interface. However, the 
fact that this interface shear stress can have an appreciable effect on 
the condensation process has been shown in reference 2. When condensation 
is taking place under the condition of a forced flow and a reduced body 
force, the interface shear stress may become the controlling force of the 
condensation process and a better understanding of gas -liquid interaction 
phenomena becomes important . 

The purpose of the present paper is to study: (l) the interaction 

phenomena between the gas -phase boundary layer and the layer of the con- 
densate which control the condensing rate and the heat transfer, and 
(2) the added effect of the body force (g force) on the condensing rate 
and the heat transfer. For the stated purpose, we shall choose a simple 
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shape, the flat plate, and. shall analyze the film condensation under the 
influence of the interface shear stress and the hody force. 

This report is divided into three parts . Part I contains a 
description of the problem and the method of analysis. Part II supplies 
the mathematical details. Part III contains a discussion and summary 
of results . It is suggested that casual readers may omit part II . 


LIST OP SYMBOLS 


D 

F(0) 

F"(0) 

f 

Fr 

Fr x 

g 

h° 

k 

L 

m 

Nu x 

Pr 

%■ 

Re 

Re x 


specific heat at constant pressure 
coefficient of binary diffusion 

Blasius stream function for gas phase at the interface, defined 
by equation (45) 

interface shear stress defined by equation (48) 

dimensionless stream function for liquid film, defined by equation 

(26) 

2 

Froude number based on L, — - — 

gL 

U M 2 

Froude number based on x, 

gx 

constant body force per unit mass 
heat of condensation 
thermal conductivity 
reference length 

mass fraction of condensable vapor in the gas -vapor mixture 


Nusselt number based on x. 


q w x 


Prandtl number. 




(T 0 -T w )k 


heat -transfer rate to plate per unit area 
Reynolds number based on L, — — 


11-y 0^ 

Reynolds number based on x, — 1 — 


v 
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Schmidt number, - 
D 

absolute temperature 

x component of velocity 
u 

^V ,00 

y component of velocity 
v 

U-V,oo 

distance along the interface 

x 

L 

distance normal to interface (see fig. l) 

y 

L 

1 c p(Tq-Tv) 

Pr h ° 

modified thickness of liquid film, defined by equation (23) 
liquid -film thickness 

similarity variable for liquid film, defined by equation (22) 

Blasius similarity variable for gas boundary layer, defined by 
equat ion ( ^3 ) 

T - T v 

To “ T w 

dynamic viscosity 
kinematic viscosity 

streamwise variable defined by equation (2l) 
dens ity 

interface shear stress 

stream function for liquid film, defined by equations (2*0 and (25) 



Subscripts 


V 

w 

o 


0 , 1,2 


00 


10 


condition of gas phase 
condition at the plate 
condition at gas -liquid interface 

order of perturbation defined by equations (37) through (39) 
free gas stream 

y Nux 

value of the parameter which would vary — — J Re„ or by 10 

u v,°o v x s/ltex 

percent from those for g = 0 


Superscript 

1 total differentiation with respect to the variable concerned 

Symbols without subscripts refer to liquid unless otherwise specified. 


I. GENERAL PROBLEM AND METHOD OF SOLUTION 
Description of Problem 


We consider that a mixture of a condensable vapor and a noncondensable 
gas is flowing over a flat plate whose uniform temperature is below the 
condensation temperature of the vapor. We further consider that the 
condensate wets the surface and that the entire system is subjected to 
an accelerative body force in the direction of flow and parallel to the 
plate. When a steady state is achieved, the following phenomena will be 
observed on the plate . A layer of condensed liquid will be flowing along 
the surface of the flat plate, dragged by the gas -phase boundary layer 
over it and also influenced by the body force (see fig. l). Condensation 
will continuously take place at the gas -liquid interface, and the heat 
of condensation will be transferred to the plate across the liquid layer. 
The present problem is then to predict theoretically the condensing rate, 
the heat transfer, and the liquid film thickness for the steady state. 

In order to make the analysis tractable, it is assumed that the flow 
of the gas and the liquid is stable and laminar, and that the condensate 
has constant properties . 
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Method of Analysis 


The thickness of the liquid layer is generally very small and, in the 
present analysis, it is considered to he in the order of a houndary- 
layer thickness . There are two layers of fluids which must he studied 
simultaneously. They are the liquid layer (liquid film) and the gas- 
phase boundary layer over the liquid. Solutions are first obtained for 
the velocity and the temperature distributions through the liquid film 
in terms of the yet unknown temperature 1 and the shear stress at the 
gas -liquid interface. The interface temperature, shear stress, and mass 
transfer are then respectively matched between the solution of the liquid 
film and the existing solution of the gas -phase boundary layer to obtain 
the final condensing rate and heat transfer. 

Liquid layer . - Let us first consider the liquid film. The governing 
equations for the liquid film are the usual laminar boundary -layer 
equations with a body force term. The equations to be solved are essen- 
tially the same as the equations used in reference 3 to study the effect 
of vehicle deceleration on melting surfaces. The equations are also 
similar to the equations used in reference k to analyze the effect of 
the buoyancy force in a boundary -layer flow. The basic differences 
between the analyses of the references and the present analysis lie, as 
far as the solution of the liquid film is concerned, in the boundary 
conditions and, particularly, in the positions at which they must be 
applied. In the analyses of the references, the boundary conditions are 
applied at the two predefined boundaries which are at finite boundary 
and at infinity. In the present analysis, the boundary conditions must 
be applied at the gas -liquid interface and at a position which is a 
liquid -film thickness away from the interface . The thickness of the film 
is not a predefined quantity and it may be obtained only as a result of 
the solutions. The thickness, moreover, varies along the plate. 

In the present study, the solutions of the governing equations are 
obtained by expanding the dependent variables and the modified film 
thickness in terms of the reciprocal of Froude number, l/Fr x . The 
solutions thus obtained with the help of a digital computer show the 
effect of the body force an a perturbation to the solution with no body 
force when a finite number of the terms of the series is retained. The 
solutions presented here are, therefore, applicable only when l/Fr^^l- 
A solution could be similarly derived which would show the effect of the 
interface shear stress as a perturbation to the solution with no interface 
drag. The present study, however, is mainly concerned with the relatively 
small body forces which exist in space flight and, therefore, the former 
solution only was considered. 


^-The interface temperature is unknown only when the gas phase is a 
mixture of a condensable vapor and a noncondensable gas. 
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The major parameters Nux/^Rex, (v o /u v ^ 00 ) N /Rex, an( l (5/x) N /Re x are 
derived from the solutions of the governing equations as functions of the 
parameters J pp F"(0), (l/Pr ) [cp(T 0 -T w )/h° ] # Froude number, and 
Prandtl number. Then for the given values of Prandtl and Froude numbers 
and the parameter, ( l/Pr ) [cp(T 0 -T w )/h° ] , the proper values of Nu x /\/Re x , 
(v 0 /u v,oo)n/R®x.> and - (8/x\/Re x for a particular problem are obtainable as 
soon as F"(0) is found. The interface shear stress represented by 
F n (0) depends on the momentum, mass, and energy interactions between 
the gas -phase boundary layer and the liquid film. The gas -phase boundary 
layer and the interaction phenomena will be discussed in the following 
section and the method of determining F"(0) will be explained. 

Gas -phase boundary layer and interface interaction .- Consider now 
the gas -phase boundary layer over the liquid film. The velocity at the 
interface is much smaller than the free -stream gas velocity, but not 
zero. In the analysis of the gas -phase boundary layer, however, the 
velocity at the interface is assumed to be zero in order that the results 
of the existing solutions of the boundary -layer equations may be utilized. 
There is a net mass transfer at the interface toward the liquid film 
because of the condensation process. The flow of gas near the interface, 
therefore, is essentially a boundary -layer flow with suction. It is 
known that suction in a boundary -layer flow increases the drag drastically 
The increased interface shear stress due to the condensation process 
itself will be shown to increase the condensing rate. The increased 
condensing rate, in turn, means again increased interface shear stress. 

It is seen then that the gas -phase boundary layer and the liquid film 
complement each other in increasing the interface shear stress represented 
by F"(0). Eventually, an equilibrium value of F"(0) will be reached 
for which both the gas -phase and liquid -phase requirements are satisfied. 
It is obvious, however, that the equilibrium value of the interface stress 
thus reached could be considerably greater than that with no mass transfer 
The equilibrium value of F"(0) is obtained in the following manner. 

The ratio of the suction to the stress functions, F(0 )/f"( 0), is first 
obtained as a function of F"(o) for the gas -phase boundary layer from 
the existing solutions of the boundary -layer equations such as are found 
in references 5 and 6. Next, from the solutions of the governing equa- 
tions, the same ratio, F(0)/F”(0), for the liquid film is derived also 
as a function of F"(0). The two functions F(o)/f"( 0) for the gas phase 
and liquid phase are then compared and the matching values of F(o) and 
F"(0) are obtained. 


II. DETAILS OF ANALYSIS 


The outline of the method of analysis has been given in part I of 
the report. In the present part, the mathematical details will be 
given . 
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Governing Equations 


Under the assumptions of the thin, stable, and laminar liquid film 
with constant properties and low speeds, the following equations describe 
the behavior of the liquid film. The equations are essentially the 
normalized boundary -layer equations, derivations of which are found in 
reference 7 • 


Continuity 


Momentum 


Energy 


where 


Su , bv _ 

5x Sy 


0 


TT 5u Tr Su _ i ^ l 

Sx bi Re ^Y 2 Fr 

U + v — = 1 

bx bi RePr Sy 2 


U = -2L_ , V = — 7_ 
u V,oo u V,oo 

e 

To - Tw 



( 1 ) 

( 2 ) 

(3) 

(M 

(5) 

( 6 ) 


The boundary conditions are as follows : 
For Y =0 (at the interface) 


3u = 

— t 0 (x) 

(7) 

cfY 

’ u -V,ooM' 

b§_ _ 

pu v ,ooVLh° 

(8) 

Sy 

(T 0 -T w )k 

e = 

1 

(9) 


8 


( 10 ) 
( 11 ) 
( 12 ) 


For Y = S/l (at the plate) 


U = 0 
V = 0 
0=0 

Boundary condition (7) is a statement of the shear stress at the 
interface, and will be evaluated by matching with that of the gas phase. 
Boundary condition (8) is derived from the consideration that the heed 
released at the interface by the condensing vapor is carried away into 
the film by conduction. Boundary condition (8) is strictly true only 
when the gas -phase boundary layer is isothermal. For practical purposes, 
however, it applies to the nonisothermal vapors also and this point will 
be explained following the gas -phase boundary -layer equations. 


For the incompressible 2 gas -phase boundary layer, we have the 
following well-known equations: 

Cont inuity 
Momentum 


Diffusion 


The boundary conditions are as follows : 

For Y = 0 

U-y = 0 (assumed approximation) (16) 

V v = g- V (17) 

v 

m = iHq (saturation concentration (l8) 

corresponding to T Q ) 


8Uy + 8Vy _ q 

8x 8 y 


Uy + v v ^ = -i 

cQC cY Re v cY 2 


2 . 

XJ + y 8m _ 1 8 m 

V 8X V 8y ScyRe v 8Y 2 


(13) 

UM 

(15) 


2 An analysis very similar to that presented here can be obtained 
for the compressible boundary layer if it can be assumed that Uy/uy ^ 

a Py,oo / Py ‘ 
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For Y = -oo 


Uv = 1 

( 19 ) 

m = ifl-oo 

( 20 ) 


When the vapor is superheated a temperature gradient exists across the 
"boundary layer, and the gas -phase energy equation is also needed to 
analyze the energy transfer associated with the superheat. The energy 
associated with the superheat, however, is very small compared to the 
latent heat of condensation for most cases. Superheat, therefore, is 
neglected in the present analysis of heat transfer and thus the energy 
equation for the gas phase is not needed. 


Solution of Equations 


Our main concern is the solution of the governing equations for the 
liquid film since, as will "be seen later, the solution of the equations 
for the gas -phase boundary layer can be obtained from the literature. 

We shall, therefore, begin the analysis with equations (l), (2), and ( 3 ). 

Experience in boundary- layer analysis tells us that a similarity 
solution of equations (l), (2), and ( 3 ) could be obtained when l/Fr x 
is either zero or infinity. The present study is mainly concerned with 
the problem when l/Fr x is rather small. We shall, therefore, handle 
the effect of l/Fr x in the form of a perturbation to the similarity 
solution for l/Fr x = 0. It was stated in part I that one of the com- 
plexities facing us in solving the governing equations is that one set 
of the boundary conditions must be applied along 6(x) which is not 
known a priori . In order to handle this complexity and also to arrive 
at the particular series form of the solutions sought, the original 
equations (l), (2), and ( 3 ) and the boundary conditions ( 7 ) through (12) 
are first transformed to more convenient forms as follows: 

We define a set of dimensionless variables as 


(21) 


■n 

Equation (22) defines the variable q in such a way that the boundary 
conditions on all the subsequent transformed equations may be applied at 
q = 0 and q = 1 only. It can be seen from equations (2l) and (22) that 
the variable A(|) is related to the liquid-film thickness S(x) as 


= _JL_ = 


5(x) A(|) 


( 22 ) 
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(23) 

The continuity equation (l) is satisfied in the usual manner by defining 
a stream function T(X,Y) by the equations 


Now we let 


*(X,Y) = A(l)f(l,Tl) 


Then equations (24-), (25), and (26) yield 




/— i— (fA ! + A — 
I ReFr \ Si 


The continuity and the conservation equations (l), (2), (3) are now 
transformed to the following set of equations: 


A 2 (S*L ^ 2f _ SLL - - AA'f — = — 

\St] Si Sri Si S n 2 J Sri 2 Sr) 3 


a 2 ^ - A (f A’ + A^]^=i-^ 


St] Si 


Si/ Sri Sr) 2 


The boundary conditions ( 7 ) through (12) become: 
For ri = 0 


U-V .00^ N RS 


Tq(X) 


(31) 
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and for t] = 1 


= lA iFr pvh° 

*1 Re k(T 0 -T w ) 

0=1 



f = 0 , 

0=0 



(32) 

(33) 

( 3*0 

(35) 

(36) 


The above transformed equations are solved by first expanding the 
variables f, 0, and A in terms of | as 


f(l,Tj) = f o ( i") ) + ifi(ri) + i^sU) + • • • (37) 

0(S,t)) = 0 o (r!) + !(ti) + i 2 0 2 U) + • • • (38) 

and 

AU) = | l/2 Ao + | 3/2 A x + g 5/£ A 2 + . . . (39) 


When the series (37), ( 38 ), and (39) are incorporated into equations 
( 29 ) and ( 30 ) when the terms with the equal powers of | are 
collected, there result the equations 


Aq 

f m , f 1 -p " = o 

I o + 2 IoI ° U 
£ d o" + foSo’ = 0 


(40) 


f 1 " - A q 2 (f 0 'f x ’ - l f 1 f 0 " - l f !% - l) + 2A 0 A 1 f 0 f 0 n = 0, 


Pr 


1 fl i" + Aq 2 (i f o 0 x ' - f o ‘0 x + | f x 0o T j + 2A o A x f o 0 o ' = 0 ' 


( 41 ) 
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Now consider the expansion of the boundary conditions. First, let us 
consider the boundary conditions (31) and (32). Before we can expand 
these boundary conditions , we must know something about the functions 
t 0 (X) and V D (x). For this purpose, we shall briefly investigate the 
nature of the solution of the gas -phase boundary -layer equations (13) 
and ( l4 ) . 

The continuity and the momentum equations (13) and (l4) are 
transformed to the following Blasius equations (see ref. 7) with the 
assumption that the dimensionless stream functions F depend on the 
variable tj^ only; that is, we have a similarity solution: 



F m (t] b ) + |F(ri B )F ,, (n B ) = 0 

(42) 

where 




/ U V , 00 

^ = v v x 

(43) 

and 

y, <v = C4 

(44) 

The boundary 

conditions (l6), (17), and (19) become: 


For ■q-g = 0 




F = ~ 2v v, 0 A-V 

V U Y jccry 

(45) 


0 

11 

h. 

(46) 

and for tj^ = 

•00 



1 — 1 
ii 

(47) 

It is seen from equation (42) and the boundary conditions (4-5) through 
(47) that F indeed depends only on rj-g if V Vj0 is proportional to 
1 /n/x and therefore F(0) is constant. In reference 5 the solution of 
equation (42) is obtained for many different values of F(0). 3 Also an 


3 There exist the following relationships between the present function 
F and the corresponding stream function used in reference 5* 

F( 0) = ~[F(0) of ref. 5] 


F"(0) = - i [F”( 0) of ref. 5] 
4 
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approximate solution of the equation can he obtained from reference 6. 

The interface shear stress r 0 (x) can be derived from the solutions as 

t 0 (x) = - (W) 

From equations (4-5) and (48) we see that when the condensing rate and 
t d (x) vary as l/x 1/2 we have similarity in the gas -phase boundary layer 
and the solution can be obtained from reference 5* 

Now let us examine the liquid film to determine whether the condensing 
rate and the interface shear stress vary with 1/x and therefore 
whether the solution of the liquid film is compatible with that of the 
gas phase. When g = 0 (g = 0), an examination of the appropriate govern- 
ing equations (4o) and the boundary conditions (31) through (36) shows 
that the condensing rate V Q and the shear stress along T] = 0 vary with 
1/sE and, therefore, they can be made to match with the solution of 
equation (42). When g > 0, on the other hand, the distribution of the 
condensing rate and the shear stress deviate from that corresponding to 
1 /n/x. However, for the small values of g considered in the present 
analysis, the deviation will be accordingly small and the effect of the 
deviation on the final results for g > 0 is neglected. The boundary 
conditions (31) and (32) are now expanded to the following form. 


f o n (0) + gf x n (0) + g 2 f 2 n (0) + 

e 0 '(°) + se i’(o) + . . . 


. = F"(o)jSte (a 0 + la, + £%, + . . .) 

(**9) 


Pr h 


c p(4'o ” T w ) l 2 


4 Ao 2 f o (0) + g 


| Ao 2 f x (0) + 2 A o A 1 f o (0) 



(50) 


Expansion of the rest of the boundary condition is rather straightforward 
and the boundary conditions (31) through (36) become for equations ( 4 o) 
and ( 4 l), respectively: 


For i) = 0 


r n _ /PvPv «n A 

f ° -J-p (r FA ° 

Aq fo _ _1_ gp(T 0 ~ T v ) 
2 (V ~ Pr h ° 

0 O = 4 


) ( 51 ) 


J 



fi* 


\ 4ofi + SAoAifo L Cp(T 0 - T v ) 


and for T) = 1 


Pr 


0i = 0 


fo = V = 0 


Qn ~ 0 


> (52) 


( 53 ) 


fi = fi* =o 
0 ! = 0 


( 5*0 


A study of equations (4o) and (4l) and boundary conditions (5l) through 
(5*0 shows that the solutions of the equations depend on the three 
parameters, Pr, J P-yPy/ PM- F"(o), and ( l/Pr ) [cp(T 0 -Ty)/h°] . Solutions of 
the equations are obtained by the use of an IBM 70A digital computer for 
several combinations of the parameters. Each system of equations, 
equations (^0) and (^l), contains three unknowns whereas it has only two 
equations. The unknowns are the functions f, 0, and the constant A. 

We have, however, one extra boundary condition for each of the systems. 

In order to solve the systems, the constant A is first assumed and the 
system is integrated. The compatibility of the result is then checked 
for the extra boundary condition. The foregoing procedure is repeated 
until the iterative scheme on A produces a solution which satisfies 
all of the six boundary conditions. A typical result is shown in figures J 
and 3- Figures 2 and 3 are plots of the solutions of equations (4o) and 
(4l), respectively, for the liquid Prandtl number of unity, 

J PyPy/ pp F"(0) of -0.01, and (l/Pr)[cp(T 0 -T w )/h°] of 0.05 . Figure 2 shows 
the velocity and temperature profiles for the case of zero body force. 
Figure 3 shows the first-order effect of the body force on the profiles. 


Heat Transfer, Condensing Rate, and Film Thickness 


For a given set of the parameters, Pr, ^ PP F n (0), (l/Pr) 
[cp(T 0 -T v )/h°], and l/Fr x , the heat transfer, condensing rate, and film 
thickness can be obtained, respectively, from the following equations. 
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_Nux 

n/R®x 




(55) 


u V,co 


sl^X 


1 

2 


A o£o + ( A cfl + A l f o) + 
- X 


(56) 


| *i^ e x - A o + jj£~ A 1 + 


(57) 


X 


Some of the typical results are shown in figures 4 through 6. Figur e 4 
shows th e dim ensionless heat transfer and condensing rate, Kux/jRe x and 
v 0 /u v oo s/Rsx? as functions of the interface shear-stress parameter 
/pipy pp F rt (0) for the case of zero "body force. There are two sets of 
curves in the figure. They are for ( l/Pr )[cp(T 0 -T w )/h°] of 0.005 and 
0.05, respectively. Figure 5 shows the dimensionless film thickness as 
a function of the same variables as figure 4. Figure 6 shows, for the 
two values of ( l/Pr ) [c- p (T 0 -T y)/h 0 ] , the value of l/Pr x which would 
cause Nux/s/Rev and v 0 /uy >0 o jRe x to deviate by 10 percent from those for 
the case of 1/Fr x = 0. Figure 6, therefore, shows the effect of the 
body force on the heat transfer and the condensing rate as a function of 
F”(0). The parameters .Jpyy pp F n (0) and ( l/Pr ) [c p (T 0 -T w )/h 0 ] 
now must be matched to the appropriate parameters of the gas -phase 
boundary layer in order that the final solution of the problem may be 
obtained in terms of the free -stream gas characteristics. The method of 
obtaining the equilibrium value of F”(o) will now be explained for the 
case of the gas phase comprised of a pure vapor. In this case, the 
entire resistance to heat transfer is essentially offered by the liquid 
film alone. We derive from equation (45) the expression 


F(0) 

F"(0) 


J 


•Vq 

uy, 

PyMy 

PM 


■ F"(0) 


(58) 


The left side of the equation, F(0 )/f"( 0), for the gas -phase boundary 
layer is obtained from data given in references 5 and 6 and is expressed 
as a function of F n (0) in figure 7* Next the right side of equation (58) 
is obtained from figures, such as 4 through 6, for a given value of Pr, 
(l/Pr)[cp(T 0 -T v )/h °], 4 and l/Fr x . The values thus obtained are then 
plotted in figure 7, against F ,f (0) for a given value of ^pyiy pp. The 
intersections, such as point a, shown then give the desired equilibrium 
value of F rt (0) for the particular problem. Once the equilibrium value 

4 The interface temperature T 0 is the saturation temperature of the 
vapor corresponding to the pressure . 
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of F n (0) is known, parameters such as Nu x /y/Re^ and v 0 /uy 2<Xs \l^ e x can he 
immediately obtained from figures 4 through 6. 

Next consider the case in which the gas phase is a mixture of a 
condensable vapor and a noncondensable gas. For this case, resistance to 
heat and mass transfer is offered both by the gas -phase boundary layer 
and the liquid film. It can be seen from reference 6 that, for Schmidt 
numbers sufficiently close to 1, the diffusion process through the gas- 
phase boundary layer requires that the following relationship must be 
satisfied in addition to equation ( 58 ) in the present problem. 

F( 0) _ n ^£0 ~ m Q 

F"(0) 1 - mo ^ 

The additional equation (59) can be satisfied because now the saturation 
temperature T 0 is an unknown quantity and is a f inaction of mQ. 

Derivation of the equilibrium value of F”(0) for this case is 
considerably more complicated. For this case, ( l/Pr ) [cp(T 0 -T w )A°1 is 
not a priori a known factor because T 0 is not known. One may start an 
iteration analysis by assuming a suitable value of To, and therefore 
corresponding values of mo and ( l/Pr) [cp( T 0 -T w )/h°] . A point such as 
a is then obtained in figure 7 in the same manner as for the case with 
the pure vapor. The value of F(0 )/F u ( 0) corresponding to a point a is 
next compared to the right side of equation (59)* If equation (59) is not 
satisfied, then another value of T 0 is assumed and the foregoing pro- 
cedure is repeated. The iteration pr ocess is conti nued until equation (59) 
is satisfied. Proper values of Nux/\/Re^, v 0 /u Y ,oo .jRe x , etc ., could then 
be obtained for the case of the gas phase comprised of both condensable 
vapor and noncondensable gas. Results of this nature will not be pre- 
sented, but can be obtained for any specific case by the procedure outlined. 


III. DISCUSSION AND RESULTS 
Discussion 


It was shown in part II that the heat transfer, condensing rate, and 
other phenomena depend on the four parameters, Pr, ! p.u / PP F”(0), 
(l/Pr)[cp(T 0 -T w )/h°], and l/Fr x . The analysis of reference 8 showed that, 
in the process of film condensation due to gravity alone, the solution 
of the problem did not explicitly depend on Pr except for extremely 
small values of Pr. For the present problem, the effect of the Prandtl 
number on the heat transfer and the condensing rate is shown in figure 8 
for p F"(0) = -0.001, (l/Pr)[c p (T 0 -T v )/h°] = 0.05 and l/Fr x = 0. 

It is seen that for the conditions represented in figure 8, the heat 
transfer and condensing rate do not explicitly depend on the Prandtl 
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number appreciably. Though the effect of Pr was checked here only for 
the particular values of the parameters, J p ^ / pp F n (o), 

( l/Pr ) [ c p ( T 0 -T w )/h° ] , and l/Fr x , it may be reasonable to believe, in the 
light of reference 8, that the solution of the present problem does not 
depend on Pr explicitly except for extreme values of the Prandtl number. 
The solution, therefore, mainly depends on the three parameters, 

Jp v p v / pp F"(0), (l/Pr)[c p (T 0 -T w )/h°l, and l/Fr x , and the effect of the 
Prandtl number appears through the parameter (l/Pr)[cp(T 0 -T w )/h°] only. 

Figure 9 shows the heat transfer and condensing rate as a function 
of J P v (-t y / pp for the two different values of (l/Pr)[cp(T 0 -T w )/h°] of 
0.005 and 0.05> respectively. The figure is for l/Fr x = 0 and for the 
condensation of a pure vapor. It is seen that the condensing rate 
increases slowly with J P T P V / p|T and increases rapidly with 
(l/Pr)[cp(T 0 -T v )/h ]. The fact that it increases rapidly with the param- 
eter ( l/Pr ) [ c p ( T 0 -T w )/h° J is expected from the nature of the parameter. 

The comparative insensitiveness of the condensing rate with respect to 
the parameter J p y M- v / PI 1 may he rather surprising at first glance. Fig- 
ure 4 shows that the condensing rate increases quite rapidly with 
N /p v q v / pp F"(0). Therefore, if the gas -phase boundary layer were unaf- 
fected by the condensing process, F rt (0) would be constant and the con- 
densing rate would increase with j p v p v / pp. as rapidly as shown in figure 
4. The gas -phase boundary layer, however, is influenced, substantially by 
the suction velocity at the interface which is created by the condensation 
process itself as was pointed out in part I. Assume that the density of 
the vapor is increased, which implies that the parameter J PyM-y/ PM- i s 
increased. The condensing rate, according to figure 4, should increase 
provided F"(0) remains unchanged. The increased vapor density, however, 
also tends to decrease the suction velocity of vapor at the interface 
and, therefore, tends to decrease F"(0). The net result is that the 
condensing rate increases with j PyPy PP but not nearly as rapidly as 
shown by figure 4. 

Next, consider the heat transfer. The parameter Nu x / N /Re x for 
(l/Fr)[cp(T 0 -T w )/h°].of 0.005 is seen in figure 9 to be considerably 
higher than th at for ( l/Pr ) [cp(T 0 -Ty)/h 0 ] of 0.05 in the upper part of 
the range of J p^p^/ pp considered. This, however, does not necessarily 
mean that the heat transfer increases as ( l/Pr ) [cp(T 0 -T w /h°] decreases. 
Besides the magnitude of ( l/Pr ) [cp(T o -T v )/h 0 ], the heat transfer depends 
also on the reason for the variation of the parameter (l/Pr)[ Cp(To-T w )/b°]. 
This may be seen more clea rly f rom the following example. Heat transfer 
can be obtained from Nu x / S /Re^ by the following equation 



(60) 
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Consider that, for ^ Pyfiy/ pp = 0.03, "the parameter (l/Pr ) [c- p (T 0 -T v )/h°] 
is decreased from 0.05 to 0.005* Figure 9 shows that Nu x /\/Rex is then 
increased by the factor of 1.53* Now, let us assume that the decrease in 
( l/Pr ) [cp(T 0 -T w )/h°] is caused "by the decrease in (T 0 -T w ). Then equation 
( 60 ) shows that the heat transfer is decreased by the factor 0 . 153 * 

Next, let us assume that the decrease in (l/Pr)[cp(T 0 -T w )/h°] is caused 
by an increase in h° . Then it is seen in equation (6o) that the heat 
transfer increases by the factor 1.53* Therefore the heat transfer may 
either increase or decrease with the parameter (l/Pr)[cp(T 0 -T w )/h°] 
depending on the reason for the variation of the parameter. The physical 
explanations for these phenomena can be made as follows. Figure 10 shows 
that the liquid -film thickness is reduced as the parameter 
( l/Pr ) [ c p ( T 0 -T w )/h° 1 is decreased. Now first consider the case in which 
(l/Pr)[c p (T 0 -T w )/h°3 is decreased because of a reduction in (T 0 -T w ). 

The liquid film is now thinner and offers less resistance to heat transfer 
to the wall. However, the heat -transfer potential (T c "•Tw) is also 
decreased. The decrease in the heat -transfer potential is such that the 
net effect on heat transfer is to decrease it even with the reduced film 
thickness. Next, consider the case in which the parameter, 

(l/Pr )[cp(T 0 -T-yO/h 0 ], is decreased because of an increase in h° . Now 
the only change, as far as the heat transfer is concerned, is the thinning 
of the liquid film. The heat transfer to the wall is therefore increased. 
It is also interesting to note that Nux/v/Re^ increases with ^ J 9-^^/ pp 
much faster when (l/Pr)[cp( T 0 -T w )/h°] is 0.005 than when it is 0 . 05 . 

The reason for this can be seen in figure 10. As J pyP-y/ pp is increased, 
the film thickness decreases much more rapidly when (l/Pr)[cp(T 0 -T w )/h°] 
is 0.005 than when it is 0 . 05 * 

Before leaving the case of l/Fr x = 0 it should be noted that the 
complementing interaction between the gas phase and the liquid film 
(explained in part I of the report) has a drastic effect in increasing 
the condensing rate and the heat transfer. The following numerical 
example may bring out this point a little more clearly. Consider pure 
steam condensing at atmospheric pressure. The value of the parameter 
^VPP is about 0.005* " When T 0 -T w « 90° F, the value of the parameter 
(l/Pr ) [cp(T 0 -T v )/h° 1 is about 0.05. If it is erroneously assumed that 
the condensation does not affect the gas -phase boundary layer, the 
correspondin g F n (0) is -0.33206 a ccord ing to reference 5* Figure i- then 
yields Nux/^Rex « 0.19 and v 0 /uoo w 0.009^* When the actual inter- 

action is considered, on the other hand, F n (0) is increased substantially 
because of the effect of the condensation on t he ga s -phase boundary layer. 
The resulting values of Nux/jRe x and v 0 /u v ^oo */R e x are obtained from 
figure 9 as 0.53 and 0.0268, respectively. The effect, of course, is not 
as large when (l/Pr)[cp(T 0 -T w )/h°l is 0.005. 
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The effect of the body force on the condensing rate and heat transfer 
is seen in figure 11. In the figure, the values of l/Fr x 'which would 
increase the condensing rate and heat transfer by 10 percent from the 
values for l/Fr x = 0 are plotted against the parameter ^ pyPy/ PI- 1 f° r 
(l/Pr ) [cp(T 0 -Tv)/ii°] of 0.05 and 0.005, respectively. It is seen in the 
figure that, for the values of ^ PyP^./ pp considered, the values of 
(l/Fr x ) 10 for the condensing rate at (l/Pr )[cp(T 0 -T~)/h°] °f 0.05 are one 
to two orders of magnitude greater than those at (l/Pr)[cp(T 0 -T w )/h°] of 
0.005* The values of (l/Fr x ) 10 for the Nusselt number are 1/2 to one 
order of magnitude greater when (l/Pr ) [cp(T 0 - Tw)/h°] is 0.05 than when 
it is 0.005* Thus for the higher values of (l/Pr)[cp(T 0 -T w )/h°] it 
requires a much larger body force to affect the heat transfer and the 
condensing rate a given relative amount than it does at the lower values. 
The relative effect of the body force on the condensing rate and heat 
transfer decreases as J P-yM-y/ pp is increased at the lower value of 
(l/Pr)[c p (T 0 -T v )/h°3 whereas it becomes rather insensitive to ^ p-yPy/ pp 
at the higher values of (l/Pr ) [cpCTo-T-y-)/^] . 

The typical velocity and temperature profiles for the liquid film 
are given in figures 12 and 13* It is seen that the velocity and tem- 
perature profiles are almost linear when l/Fr x = 0. As the body force is 
added, however, the velocity profiles have a parabolic form and, at the 
same time, the liquid film becomes thinner. 


Summary of Results 


The laminar film condensation on a flat plate under the influence of 
the gas -liquid interface shear and the body force (g force) has been 
analyzed . 

It is shown that the major governing parameters in the condensation 
of a pure vapor are the ratio of the product of the density and viscosity 
of the gas to that of the liquid, ( PyPy./ pp ) 1 /z , the ratio of the tem- 
perature potential to the heat of condensation, ( l/Pr ) [c p (T 0 -T w )/h°] , 
and the Froude number, Uy,oo 2 /gx. When the gas phase is a mixture of a 
condensable vapor and a noncondensable gas, the gas -phase diffusion 
potential also becomes an important parameter. 

The heat transfer, the condensing rate, and the liquid -film 
thickness are calculated for the case of a condensing pure vapor. The 
results show that the condensing rate increases slowly with the density- 
viscosity ratio, p y p y /pp, and rapidly with the ratio of temperature 
potential to heat of condensation, ( l/Pr ) [cp(T 0 -T w )/h° ] . The Nusselt 
number also increases with the density -viscosity ratio but the rate of the 
increase becomes smaller as the ratio of driving temperature potential 
to heat of condensation increases . 
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It is found that as the condensation proceeds, it creates a suction 
velocity at the interface for the gas -phase boundary layer which results 
in an increased interface shear stress . The increased shear increases 
the condensation rate by increasing the liquid velocity. The increased 
condensing rate, in turn, increases the interface shear stress imposed by 
the gas -phase houndary layer. Because of this complementary interaction 
between the gas and liquid, it is found that the film condensation in a 
boundary -layer flow occurs at a high rate. 

The ratio of temperature potential to heat of condensation is found 
to have a large influence on the magnitude of the body force which can be 
supported by the liquid film before the characteristics of the film are 
altered substantially. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 31 , 1961 
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Figure 2-- atypical solutions of the zeroth order perturbation equations 
equation (i+O), for the stream function, v elocity p rofile, and 
temperature profile; Pr = 1.0; P = 0.05; J PyP-y/ PP F n (0) = -0.01. 





Figure 3 .- Typical solutions of the first-order perturbation equations 
equation (^l), for the stream function, v elocity p rofile, and 
temperature profile; Pr = 1.0; T = 0.05; F n (0) = -0.01. 
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Figure 13 •- l?ypical tempera ture prof ile for liquid film; Pr 

r = 0.05; ^pyi^/pp f"(o) = -o.oi. 
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